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ABSTRACT 

We present a model for steady state winds of systems with super-Eddington lumi- 
nosities. These radiatively driven winds are expected to be optically thick and clumpy 
as they arise from an instability driven porous atmosphere. The model is then applied 
to derive the mass loss observed in bright classical novae. The main results are: 

(i) A general relation between the mass loss rate and the total luminosity in super- 
Eddington systems. 

(ii) A natural explanation to the long duration super-Eddington outflows that are 
clearly observed in at least two cases (Novae LMC 1988 #1 & FH Serpentis). 

(iii) A quantitative agreement between the observed luminosity evolution which is 
used to predict both the mass loss and temperature evolution, and their observa- 
tions. 

(iv) An agreement between the predicted average integrated mass loss of novae as a 
function of WD mass and its observations. 

(v) A natural explanation for the 'transition phase' of novae. 

(vi) Agreement with 77 Carinae which was used to double check the theory. The 
prediction for the mass shed in the star's great eruption agrees with observations 
to within the measurement error. 

Key words: Radiative transfer — hydrodynamics — instabilities — stars: atmo- 
spheres — stars: individual LMC 1988 #1, FH Ser, 77 Carinae — novae, cataclysmic 
variables 



1 INTRODUCTION 

The Eddington luminosity is the maximum luminosity 
allowed for a stationary spherical, homogeneous, non- 
relativistic and fully ionised system. If one allows motion, 
then a steady state super-Eddington configuration generally 
does not exist unless the system is just marginally super- 
Eddington or it has a very high mass loss rate. Neverthe- 
less, nature does find a way to construct steady state con- 



atmosphere necessarily has to reside at significantly deeper 
optical depths. The inconsistency arises because the sonic 
point and the critical point have to coincide in a steady 
state solution, while nowhere within the plausible space of 
parameters such a solution exists. 

A solution was proposed in which the atmosph ere of 
rj Carinae is inhomogeneous, or porous ( Shaviv, 2000| ). The 
inhomogeneity is a natural result of the instabilities of atmo- 



figurati Dns in which super-Eddington luminosities exist with 
only a relatively small mass loss rate, i liis was perhaps best 
demonstrated with 77 Carinae's giant eruption. 

ri Carinae was clearly super-Eddington duri ng its 20 



year long eruption (s ee for instance the review by Davidson 
&: Humphreys, loo's). Yet, it was shown that the observed 



mass loss and velocity a re inconsisten t with a homogeneous 
solution for the wind ( ghaviv, 2000 ). Basically, the sonic 
point obtained from the observed conditions necessarily has 
to reside too high in the atmosphere, at an optical depth of 
only ~ 1 to ~ 300, while the critical point in a homogeneous 



spher es that are close to the Eddington luminosity ( shaviv 
2001). The inhomogeneities, or 'porosity', reduce the effec- 
t ive opacity an d increase the effective Eddington luminosity 
( ^haviv, 199^ ). 

In this paper, we are interested in understanding the 
wind generated in cases in which the luminosity is super- 
Eddington. To do so, it is advantageous to find a class of 
objects for which better data than for 77 Car exists. One 
such class of objects is novae. 

Novae are a very good population to analyze in order to 
understand the behaviour of steady state super-Eddington 
winds. The main reasons are: 
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(i) The mass and luminosity are better known than for 
many other objects. For example, although rj Carinae was 
clearly super Eddington, it is not clear by how much it was 
so: Its mass can be anywhere between 60Mq and IOOA/q 
and its luminosity during the eruption is even less accurately 
known. The ejected mass could have been between IMq and 
3Mq. This is not accurate enough for our purposes here. 

(ii) The opacity where the sonic point is expected to be lo- 
cated is governed by Thomson scattering. In AGB and post- 
AGB stars that generate strong radiatively driven winds, the 
opacity is a very sensitive function of the local parameters of 
the gas at the sonic point. Thus, even though their observed 
mass and luminosity can be fairly accurately deduced from 
the observations, the modified Eddington parameter which 
should take into account the opacity of the dust, for exam- 
ple, is not known reasonably well. 

(iii) Since the luminosity during the super-Eddington 
episode of novae eruptions can be significantly above the 
Eddington limit, the inaccuracy of F = L/£,Edd is less criti- 
cal to the exact value of the mass loss. In objects that shine 
very close to the Eddington limit for a long duration, with a 
relatively small mass loss rate, the theory cannot give pre- 
cise predictions. Thus, if for example the winds of the most 
luminous WR stars arises because the objects are marginally 
super-Eddington, it would be hard to compare their obser- 
vations to the theory presented here because the accuracy 
in the determination of F — 1 will be rather poor. 

(iv) Novae generally exhibit a 'bolometric plateau' in 
which the bolometric luminosity decreases slowly over a rela- 
tively long duration, significantly longer than the dynamical 
time scale. Therefore, if this luminosity is super-Eddington, 
then clearly a steady state model for the super-Eddington 
flow should be constructed. This is clearly the case with two 
specific novae: Nova LMC 1988 #1 and Nova FH Serpentis. 
We exclude from the discussion here the very fast novae for 
which this property is least pronounced. 

The fact that the 'bolometric plateau' is sometimes ob- 
served to be super-Eddington does not presently have any 
good theoretical explanation. The steady state burning of 
the post maximum of novae is often predicted to be given 
or at l east approximate d by the core-mass luminosity re- 
lation (Paczynski, 197C). The classical core-mass luminos- 



ity relation increases monotonically with the mass of the 
WD and saturates at the Eddington hmit. It does not yield 
super-Eddington luminosities. The observations of super- 
Eddington luminosities over durations much longer than the 
dynamic time scale, lead us to the hypothesis that bright no- 
vae must have a porous atmosphere. A porous atmosphere 
with a reduced effective opacity will naturally give a core- 
mass luminosity relation that increases monotonically with 
mass beyond the classical Eddington limit, providing the 
arena for steady state super-Eddington winds. 

It is these winds which are the subject of this paper. In 
section ^ we present the 'wind theory' for super-Eddington 
atmospheres. In section ^, we apply the wind theory to two 
specific novae that were clearly super-Eddington over a long 
duration, to the nova population in general and to Carinae 
and compare the results with observations. 



2 THE FUNDAMENTAL STRUCTURE OF 
SUPER-EDDINGTON WINDS (SEWS) 

2.1 Some general considerations 

What have we learned from Carinae? i] Car has shown 
us that an object can be super-Eddington for a duration 
much longer than the dynamical time scale while driving a 
wind which is significantly thinner than one should expect 
in a homogeneous wind solution. When one tries to con- 
struct a steady state wind solution, one has to place the 
sonic point at the critical point — where the net forces on 
a mass element vanish. If a system is in steady state and 
super-Eddington, then the critical point has to reside where 
alternative means of transporting the energy flux, namely by 
convection or advection with the flow, become inefficient. 
This point however, happens relatively deep in the atmo- 
sphere (r S> 1), implying that the mass loss M — A-kFI? pVg 
(where Va is the sound velocity), is very large. In f act, the 
mass loss rate becomes of order [L — C-E^dA) / (e.g., 3wocki 
h Gayley, 1997, Shaviv, 2000). However, if the radius of 



the system is fixed, then because a minimum energy flux 
of MGMi, / R = Mvcsc/'i (with Vesc being the escape veloc- 
ity) has to be supplied in order to pump the material out 
of the gravitational well, one obtains that unless Vs ^ Ucsc, 
t he radiation will not b e able to provide the work needed 
( lOwocki fc Gayley, 1997| ). This implies that below a certain 
luminosity and radius of the system, there is no steady state 
solution. Clearly, the system would try to expand its outer 
layers (that are driven outward but cannot reach infinity), 
thereby reducing the escape velocity, until a steady state can 
be reached. 

One would expect that as time progresses, the sonic 
point of the wind would move monotonically downwards, 
pushing more and more material upwards thereby expand- 
ing the atmosphere and accelerating more mass until all the 
available luminosity would be used-up to pump material out 
of the well. If this expectation is realized, r) Car would have 
appeared as a faint object with high mass loss at low ve- 
locities. In other words, observations of rj Car suggest that 
down to the optical depth of at most ^ 300 , which is th e 
deepest that the sonic point could be located (Shaviv, 2000), 
the total mass is ^ 0.02Mq, implying that this part of the 
envelope should have been in steady state for time scales 
longer than ~ 3 months. This hypothetical steady state is 
inconsistent with r) Car trying to reach an equilibrium in 
which most of the radiation is used up to accelerate a very 
large amount of mass to low velocities, the star has been 
notably super-Eddington for a long duration without accel- 
erating more and more mass at lower velocities. 

We will soon show that novae, at least during their 
'bolometric plateau' defy the Eddington limit. In some cases 
at least, a steady state super-Eddington configuration is 
reached in which the kinetic energy in the flow plus the rate 
in which gravitational energy is pumped is only a small or 
moderate fraction of the total radiative energy flux at the 
base of the wind. This too is inconsistent with a sonic point 
located deep inside the atmosphere. 

So, how did T) Car ci rcumvent its bloating up? It was 
proposed by Shaviv (200C) that the solution to the problem 
is in having a porous atmosphere. In such an atmosphere, 
density perturbations naturally reduce the effective opac- 
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ity «eff (Shaviv, 199J). Consequently, the radiative force is 



decreased and the effective Eddington limit is increased to: 



-Co 



Kcff 



> £e 



(1) 



where acq is the microscopic opacity. When kq = KThomson, 
the Eddington limit £Edd corresponds to the classical Ed- 
dington limit. In most other cases, where the microscopic 
opacity kq is larger than the Thomson opacity, -CEdd cor- 
responds to a lower modified Eddington luminosity. In the 
rest of this work, we shall not explicitly state whether the 
Eddington limit is 'classical' or 'modified'. The distinction 
should be made according to the underlying microscopic 
opacity, which is close to that of Thomson scattering in the 
specific cases studied here. We shall however make the im- 
portant distinction between the 'modified' and the 'effective' 
Eddington luminosity, the latter being the effective Edding- 
ton luminosity in a non-homogeneous system. 

A mechanism which converts the homogeneous lay- 



ers into inhomogeneous was suggested by Shaviv (2001) 



It is demonstrated that as the radiative flux through the 
atmospheres surpasses a critical Eddington parameter of 
Tcrit ~ 0.5 — 0.85 (the exact numerical value depends on 
the boundary conditions), the atmosphere becomes unsta- 
ble to at least two different instabilities, both of which op- 
erate on the dynamical time scale (namely, the sound cross- 
ing time of a scale height). Consequently, as the radiative 
ffux approaches the Eddington limit and surpasses it, the 
radiation triggers the transition of the atmosphere from a 
homogeneous one to an inhomogeneous one. The inhomo- 
geneities increase the effective Eddington limit thereby func- 
tionally keeping the radiation flux at a sub-Eddington level; 
namely, even if F > 1, the optically thick regions experience 
a Feff < 1. Other instabilities that operate in more complex 
environments could too be important and co ntribute to this 



transiti on. F or example, s-mode instabilities (Glatzel & Kiri 
akidis, ),993|) or the instability of dynamically detached outer 
layers ([Stothers fc Chin, 1993) operate under more complex 



opacity laws. The instability of 'photon bubbles' app ears 
when strong magnetic flelds are present ( Arons, 1992| ). In 



principle, since the adiabatic index approaches the critical 
value of 4/3 for instability as the Eddington limit is ap- 
proached, many mechanisms which are otherwise unimpor- 
tant do become important. 

The 'porosity' reduces the effective opacity only as long 
as the perturbations are optically thick. Therefore, a wind 
will necessarily be generated from the regions in which the 
perturbations become optically thin, since from these re- 
gions upwards the effective opacity will be the normal 'mi- 
croscopic' one and the effective Eddington limit will return 
to be the classical one. 

The Fcrit < F < 1 case: — The details of the geometry 
(or inhomogeneities) of the regions depend on the instability 
at play, and for example can be in the form of 'chimneys' 
or 'photon bubbles'. The lowered opacity is achieved by fun- 
neling the radiation through regions with a much lower than 
average density. These regions can be super-Eddington even 
when the mean F parameter is smaller than unity. In such 
cases, mass loss should be driven in the super-Eddington 
'chimneys' of lower density. These atmospheres are then ex- 
pected to be very dynamic. Once any accelerated mass el- 
ement reaches the optically thin part of the atmosphere, it 



will start to experience an average force that is at a sub- 
Eddington value and the mass flow will stagnate, proba- 
bly forming something which looks like 'geysers'. (It is very 
unlikely that the escape velocity will be attainable in the 
'chimneys' since shocks would probably limit the flow to ve- 
locities not much larger than the speed of sound). A wind 
could then be generated in the optically thin part of the 
atmosphere through the standard li ne driving mechanism 



Pauldrach et al., 1986 



( pastor et al., 1975 , or for example 
and references therein), with the notable consideration that 
the base of the wind is clumpy. 

The F > 1 case: — When F > 1 on the other hand, 
a continuum driven wind has to be generated. The reason 
is clear. Since the perturbations have to be optically thick 
to affect (and reduce) the opacity, at a low enough optical 
depth one would expect to return back to a super-Eddington 
ffow. What is this optical depth? If we climb up the thick 
yet porous atmosphere, since it is effectively sub-Eddington, 
the average density will decay exponentially with height. At 
some point, the density will be low enough that a typical 
perturbation 'element' becomes optically thin. Since the per- 
t urbations are expected to be of order a scale height in size 
( ^haviv, 2001 ), this depth would be where a scale height has 
an optical width of order unityQ. Beyond this point, the typ- 
ical perturbation on scales of order the a scale height cannot 
reduce the opacity. This is the place where the sonic point 
should be located for a steady wind to exist. 

More specifically, if the flux corresponds to an Edding- 
ton parameter F, then the optical depth at which pertur- 
bations cannot reduce the effective Eddington parameter to 
unity should scale with F — 1. The reason is that the de- 
crease of the luminosity by the effective opacity should be 
proportional to the deviation of the actual luminosity from 
the Eddington one. Namely, when close to the Eddington 
limit, a blob with the same geometry needs a smaller den- 
sity fluctuation and with it a smaller change in the opacity, 
to reduce the effective opacity by the amount needed to be- 
come Eddington. Thus, when closer to the Eddington limit, 
the sonic point can sit higher in the atmosphere. We shall 
elaborate this point to show that this is indeed the case in 
^3 



2.2 The Structure of a Steady SEW 

The above considerations lead us to propose the following 
structure of a super-Eddington wind (hereafter SEW). Con- 
sider Fig. |l| for the proposed structure of a super-Eddington 
atmosphere (one with Ltot > -CEdd) and the SEW that it 
generates. Four main regions can be identifled in the atmo- 
sphere and they are: 

• Region A: A Convective envelope - where the density 
is sufficiently high such that the excess flux above the Ed- 
dington luminosity is advected using convection. The radia- 



Note that if the atmosphere is static and therefore has an expo- 
nential density profile, this location would also correspond to the 
place at which the optical depth is of order unity. Since a thick 
wind is expected to form, the physical depth where the optical 
width of a scale height is of order unity does not correspond to 
the physical location (the photosphere) where the optical depth 
is of order unity. 
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Figure 1. The proposed structure of a super- Eddington atmo- 
sphere (one with Ltot > •CEdd) ^^'^ the wind that it generates. 
The four regions, which are described in detail in the text, are: 
Region (A) A Convective envelope in which the radiation is sub- 
Eddington - L^ad < -^Edd < ^tot- Region (B) A porous atmo- 
sphere in which the effective Eddington luminosity is larger than 
the classical Eddington luminosity: £Edd < ^rad = itot < ■Ccft- 
Region (C) of an optically thick, continuum driven wind, where 
perturbations are optically thin and the effective Eddington limit 
tends to the classical value. Region (D) of a photosphere and 
above. 



tive luminosity left is below the classical Eddington limit: 
^rad <_^Edd < itot ■ It was shown by Joss, Salpeter & Os- 



triker (1973) that convection is always excited before the Ed- 
dington limit is reached. Thus, if the density is high enough 
and the total flux is super-Eddington, this region has to ex- 
ist. 

• Region B: A zone with lower densities, in which con- 
vection becomes inefficient. Instabilities render the atmo- 
sphere inhomogeneous, thus facilitating the transfer of flux 
without exerting as much force. The effective Eddington lu- 
minosity is larger than the classical Eddington luminosity. 
>CEdd < irad = ^tot < CcS ■ V Car has shown us that the 
existence of this region allows fo r the steady state outflow 
during its 20 year long eruption ( ^haviv, 2000 ). 

• Region C: A region in which the effect of the inho- 
mogeneities disappears and the luminosity is again super- 
Eddington. When perturbations arising from the instabili- 
ties, which are expected to be of order the scale height in 
size, become transparent, the effective opacity tends to the 
microscopic value and the effective Eddington limit tends 
to the classical value. At the transition between (B) and 
(C), the effective Eddington is equal to the total luminosity. 



This critical point is also the sonic point in a steady state 
wind. Above the transition surface, Ltot > jCcB jC^dd and 
we have a super sonic wind. This wind is expected to be 
optically thick. 

• Region D: The photosphere and above. Since the wind 
is generally thick, the transition between regions (C) and 
(D) is far above the sonic surface. 



2.3 The location of the critical point 

Paramount to the calculation of the mass loss rate in super 
Eddington systems is the location of the critical point (sep- 
arating regions (b) and (c) in Fig. |l|) , which corresponds to 
the sonic point of the wind generated at steady state. There- 
fore, we should try and estimate its location and density. Its 
definition is the location where the net average force on a gas 
element vanishes, or in other words, it is where the effective 
Eddington parameter Feff is uruty. To calculate it, we need 
to know the functional behavior of the effective opacity as a 
function of height. 

The general behavior clearly depends on the charac- 
teristics of the nonlinear state of the porous atmosphere. 
However, this state is still part of an open problem under 
investigation and is therefore unknown. It will be shown, 
nevertheless, that the unknown geometrical factors that de- 
termine the mass loss can be pin-pointed while some other 
unknowns, such as Fcff , actually cancel out to first approxi- 
mation. 

The effective opacity deep inside the atmosphere re- 
duces F to sub-Eddington values. However, higher up in the 
atmosphere where the perturbations become optically thin, 
the reduction in effective opacity diminishes to zero. In the 
thin limit, the effective opacity returns back to its micro- 
scopic value. We therefore write the opacity (per unit mass) 
of the medium as: 



Kofi (At 



Foff(Ar) 



-)^(Ar), (2) 



where At is the equivalent optical width of a pressure scale 
height if the medium had been homogeneous. Ft,ff,oo is the 
effective F in the limit of large optical depths. It is a func- 
tion of F which is unknown at this stage since we don't have 
a model for the nonlinear behavior of the inhomogeneities. 
The function Feff(Ar) should be obtained from a non lin- 
ear model (which is part of work in progress). The function 
'0(At) gives the functional reduction of the opacity as a 
function of At. It should depend on the geometrical charac- 
teristics of the inhomogeneities. For At — > 0, — * 1 and 
will approach the microscopic value of the opacity kq. For 
large optical depths, we have At — + oo, and obtain ■(/) ^ 
and Kcff (Foff,oo/F)Ko. 

Next, the sonic point for a steady-state transonic flow 
must be located at the critical point, where the local 
Feff(Ar) is unity, i.e., where (ftcff (Ar)/Ko)F — 1. We there- 
fore have: 



iA(Arsc 



(3) 



Let us assume that for small optical depths, i/)(Ar <^ 
1) ~ 1 — AAt^' and that for large optical depths, ^1){At 2> 
1) ~ BAt'^*^ . If we know the constants A and B and the 
powers pi and ph, we can estimate the average density at 
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the sonic point. To do that, we need a relation between the 
average density po and Ar. To this goal we write the pressure 
scale height Ip as 



1 



(4) 



The constant u relates the effective speed of sound to the 
adiabatic one (vs). If for example the atmosphere is isother- 
mal, we have v = I/7. If its temperature gradient is derived 
from radiation transfer and it has a constant opacity, we 
find u = 3/(47). The factor 1 — Feff^oo originates from the 
'pufBng-up' of the atmosphere due to the radiative force 
(which reduces the effective gravity). Let the average den- 
sity be po, then the optical width of a pressure scale height 
if the layer is homogeneous is given by: 



Ar — pKolp 



therefore: 



(5) 



p = PunitAr(l — Tcff.oo) with p^^it = g/i^novvi). (6) 
For the small optical width limit, where 1 — 1/; ^ 1, one 

has: 

Psonic(Ar <C; 1) ~ PunitArsonic(l — Teff.oo) 



(7) 



Punit(l — Tc 



1 F - 1 

AT — Feff.o 



However, small At's are obtained when 1 — -i/; ^ 1 or 
F — 1 <C F — Feff,oo, namely, only when F is very close to the 
Eddington limit. Under such conditions: 

p,onic(AT « 1) ^ punit^ (F - if ^ ( 1 - Feff,oo ) ' " 'Z"' • (8) 

API 

The opposite limit of large optical depths corresponds 
to t/j <C 1, and consequently to F 3> 1 as well. In this limit 
we have: 



Arsonic(AT»l) 

or 



1 — Foff,oo 1 
F — Feff.oo B 



(9) 



Psonic Punit ( 1 Tcff' ,oo ) ^"^sonic 

~ PunitZ?'/^'^ (F - Feff.oc)'/''^ (1 - V^fi^^f-^'^'^ 



- cff .ck: ) 

Pur.itB^'"'^ (F-l)i/'"-(l-F, 



(10) 



where for the last approximation we used F 2> 1 ^ Fcff.oo- 

Clearly, for the particular case of = p; = 1 we obtain 
that both in the large and small optical depth regimes: 



Psonic ~ Punit (.4 ^ Or B){V - 1) 



(11) 



It is interesting to note that for this particular set of power 
laws, both regimes do not depend on the value of Fcft,cx) 
Eis it cancels out: The factor 1 — Fcff,oo introduced by the 
puffing up of the atmosphere cancels the factor (1— Fcft,oo)~^ 
which appears because when the atmosphere is closer to the 
effective Eddington limit, smaller changes in the effective 
opacity are needed to reach the critical point. 

We now proceed to show that the above simple power 
law scalings are indeed obtained under both limits. Since 
the arguments given are somewhat heuristic, the function 
^(Ar) is also calculated analytically in the appendix for 



two simple cases under both limits, giving the same power 
laws. To strengthen the argumentation even more, we also 
use published results for an inhomogeneous system with a 
particular type of random statistical distribution (that of a 
Markovian statistics) for which the opacity (and -0) can be 
calculated rigorously. 

To find the power laws, we assume for simplicity that 
the medium is divided into regions of higher density ph = 
p(l -I- 5) and volume Vh = oVtot and regions of lower density 
pi = p(l — (5) and volume Vi = (1 — Q)Vtot. From conservation 
of mass it follows that phVh + piVi = poKot or p = po/(l + 
(2a- 1)5). 

The optically thick limit: The optically thick limit can 
be generally described by Fig. ^ The reduced effective opac- 
ity is obtained by concentrating mass in dense 'blobs' and 
forming rarefied medium in between through which most 
of the radiation flux is funneled through. Simplistically, we 
have high density regions with a volume Vh, a density ph 
and a reduced flux Fh and low density regions with Vi, pi 
and higher typical fluxes Fi. 

In the absence of inhomogeneities, we would have had 
an average density po and a flux Fq. For a given temperature 
gradient, we typically have Fh,i ~ Fopo/ph,i- The result is 
accurate if the perturbations in the medium have a large ver- 
tical to horizontal aspect ratio, such that the temperature 
perturbations are negligible. Under more general perturba- 
tion geometries, this approximation is only a rough estimate 
but it should give the typical change expected in the fluxes. 

The total average force / acting on a volume which con- 
tains many inhomogeneous elements, is / — Ko{VhphFh + 
VipiFi) « {Vi + Vh)nopQFo = /o, where /o is the force cor- 
responding to the homogeneous case. That is, the average 
force does not change for a given flxed average temperature 
gradient. But the flux through the system is 



F 
Fh 



VhFh + ViFi 



FoVtot 
-[{2a- 



Ph pi 



Po 



52 



(Ko/'Jcff) > 1- 



(12) 



We find an increased fiux without changing the driving tem- 
perature gradient. Said in other words, the same tempera- 
ture gradient drives now a higher fiux — the radiative con- 
ductivity increases or the effective opacity decreases. 

What happens though when the optical depth of the 
perturbations is finite? If we look at Fig. |^ we see that 
there exists a flnite volume at the interface between the 
high density regions and the low density ones with a to- 
tal cross section of about one optical depth. The material 
in this volume sees radiation from both the high density re- 
gions and the low density regions, so its flux will be averaged 
Javr ~ {Fh + Fi)/2. The amount of mass that witnesses this 
averaged flux is of order the surface area S of the interface 
between high and low density regions times pKo- Half of the 
mass will be on the dense side (labeled '1' in the flgure) 
and occupy a volume of order Vi « S/2phKo and the other 
half on the rarefled side ('2') and will occupy a larger vol- 
ume V2 ~ S/2piHQ. When calculated, the total force /tot for 
a given temperature gradient remains unchanged, but the 
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flux will not be as high: 
'Vh 



F 

To 



^ Vi 



{V2 



V2 , V1 + V2 

2 

Pi)Pq 



Pi 

- yi){ph - 



Po 

F'tot,A.T — >CXD 

Fo 2Vtot phPi 

'^eff ,00 (l-<52)2K0P0Kot 











\ph 





(13) 



If we further assume that F 3> 1, then the signiflcantly 
under-dense regions carry the lion share of the flux while the 
denser regions transport only a minor fraction. Under such 
conditions one can neglect 1/ph when compared with 1/pi 
or assume that 5 is close to unity. We then obtain: 



F 



S 



' Kopo4(l - 
The effective opacity is then: 

'I'cfF y Favr ^cff'.oc . S 



+ 



(14) 



(15) 



fo -Ftot Ko ' 4(1 — Q!)2/toPoKot ' 

The more elaborate treatment in the appendix gives that 
the prefactor (again, for vertically elongated perturbations) 
is not 1/(4(q - if) but 3/(32Q(a - 1)^) instead (i.e., 3/4 
times smaller for a — 1/2). 

The term (5/2^0^0 Kot) can be related to the optical 
depth of the perturbations. If the physical size of the per- 
turbations is d, then we can write S — HKot/rf, where H 
is an unknown geometrical factor. The geometrical factor is 
equal for example to 1 if the geometry is that of sheets of 
width d. It is higher for higher dimensionally of the struc- 
ture (e.g., H « 2 or 3 for cylinders or spheres of diameter d, 
ff a = 1/2). 

If we now compare this result with the definition of the 
function 1/), we have: 



3S 



1 1 

32a(l - Q)2KoPoKct " \32ail-ayj ]3A^ 



3H 



,(16) 



where /3 = d/lp. It is the ratio between the typical size of 
the perturbations and the pressure scale height of the atmo- 
sphere. Since the typical wavelengths which are the quickest 
to grow are those with a size comparable to that of the scale 
height ( ^haviv, 2001 ), this factor should too be of order 
unity but could be somewhat larger or smaller, depending 
on the details. The expression in parenthesis describes a ge- 
ometrical factor which is of order unity or somewhat larger 
(e.g., if a = 1/2 then it is about 1-3 for typical H's, and it 
would typically be larger for other a's). 

Through the definition of ip in the thick limit, we have 



B = 



3H 



32a/3(l 



and Ph = 1. 



(17) 



The optically thin limit: We assume for simplicity that 
the medium is made up of spherical regions of higher den- 
sity and those with lower density, each occupying a similar 
volume. Under this approximation, we can estimate both A 
and pi without cumbersome algebra. 

The equation of radiative transfer for a gray atmosphere 

is: 



d/(x) 

ds 



= x(x)(B(x)-/(x)), 



(18) 




Figure 2. The effects of a finite optical depth on the flux and 
opacity of a system close to the thick limit. When the perturba- 
tions are optically thick, the fluxes inside and outside of the dense 
regions do not 'mix'. However, because of the nonlocal nature of 
radiation transfer, a 'blob' with a flnite optical width wiU have 
a region around its interface with the low density gas in which 
the radiative flux is averaged. Part of this interface region (de- 
noted by '1') has the high density, and part of it (denoted by '2'), 
has the lower density. The optical width of this region is of order 
unity. The result is that a finite fraction of the system's volume 
is not effective in reducing the opacity. 



where x ~ ^^P a-nd B are the opacity per unit volume (i.e., 
the extinction) and radiation source function respectively. 

The formal solution to the radiative equation for a point 
X located many optical distances from the boundary is (e.g., 
Mihalas fc Weibel Mihalas, 1984|): 



7(x, n)= / x(x — ns)_B(x — ns)e 



/o X(x~ns')ds' 



ds. (19) 



Next, we calculate the intensities /(xA,n) at point 'A' 
located inside a spherical region of radius vq with a density 
po{l + S) (cf fig. m). We assume that outside the given high 
density region the photons seen many high and low density 
regions and hence an average density. In other words, there 
is a high density correlation inside a given blob but it is 
assumed to be negligible for distances larger than the blob. 

We also assume that the radiation source function is 
B = Bo + bz and that there is no contribution from the 
perturbations. In the appendix, it will be shown to be a valid 
assumption because the correction to the source function is 
of second order, i.e., it falls as [roXo)^- We also neglect higher 
order terms in z. We define now p as the angle between the 
integration path and z. For simplicity, we also assume that 
Q = 1/2. 

Using the formal solution, the radiation intensity is then 
given by: 

/(xA,n) = Bo+Xocos(^i) (20) 

I 

Xos(l + <5) exp(— xos(l -I- S))ds 



+ 



i exp(-xoPo((l + S)ro + (s - ro)) 



Bo + —{l~SroXo)cosp + 0{roXor- (21) 
Xo 



The flux is then given by: 



1 



Fh^FA^-i> H^A, n)ndn = Fo(l - 5) 



(22) 



where Fo — 47r6/(3cxo). 
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The flux in the under-dense regions is obtained by the 
transformation S —S. In the simplest case discussed here 
of both regions having the same volume fraction, we see 
that the total flux through the system remains unchanged. 
However, the mean force is now reduced: 



/ = VhphHoFh + VipiKoFi = /o(l - roKopoS ). 



(23) 



The effective opacity is therefore reduced as well. Using the 



result of eq. ([12 
paring to eq. (T 
get 



for the thick limit and a = 1/2, and com- 
) with iP{At < 1) ^ 1 - AAt^', we finally 



A- 



2Ar 



^ = /3/2 and pi = 1, 



(24) 



with d = 2ro. In other words, also in the thin limit, the 
function depends only on geometrical factors (and not for 
example on Feg^oo). 

A Markovian mixture: A special case which can be solved 
more rigorously for any optical thickness is that of a Marko- 
vian mixture of dense and rarefied phases. A Markovian mix- 
ture is defined through the spatial statistics of the distribu- 
tion: On any ray passing through the medium, the interface 
between material (or phase) A and material B form a Pois- 
son process. In other words, if at point s the fluid is of type 
A, then at an infinitesimally close point s -I- ds, the proba- 
bility of finding fiuid B is ds/Ao. If point s is in fluid B, then 
the equivalent probability of flnding fluid A at point s -f- ds 
is ds/Xi. The effective transport p roperties in such a mix- 



Le yermore et al, 
§3.4 of l^omraning ( 



1992 



ture were extensively analyz ed by Levermore et al. (1986) 
and [Le vermore et al. (1988 ) , and it is described in detail in 

It was found that the effective 



absorption cross-section in the absence of scattering is: 



0"cff ~ 



{Poo-Q +picri + aoaiXc) 
(1 -I- {poai +pi(7o)Xc) 



(25) 



with Ac = AoAi/(Ao -f- Ai), Pi = Xi/{Xo + Ai) and ao,i being 
the two absoption cross-sections of materials A and B re- 
spectively. To translate this result to our problem notation, 
we identify po with our q, pi with 1 — a and Aq with d, such 
that Al is d(l — a)/a and Ac = (1 — a)d. Moreover, we can 
identify cro with ph^o and cri with pikq. We then flnd that 
the effective opacity at the large optical depth limit is: 



«eff,< 



KQ 



(1 + (2q- 1)5)2 



(26) 



Next, we write d — pip and At — IppoKo- Then, if we further 
assume that this reduction is large, as we did before, we can 
assume that ph ^ Pi or 5 ~ 1. Using this result, in the two 
optical depth limits {At <^ 1 for A and At ^ 1 for B), it 
is straightforward to show that: 



^ = (1 - a)f3 and B = 



(l-a)/3' 



(27) 



as well as ph = Pi = 1. Namely, in the special case of a 
Markovian mixture, we flnd that the function i/)(Ar) be- 
haves in the same way as in the simple models presented 
before, with only the normalization constants A and B be- 
ing somewhat different but of the same order as before. 

Although the analysis carried out here for the location 
of the critical point in both limits is far from accurate, it 




Figure 3. The effects of a finite optical depth on the fiux and 
opacity of a system close to the thin limit. Point A located in 
a higher density region will witness a slightly reduced flux from 
the average one. The main reason is that when calculating the 
intensities I{xa, m) at point A, the physical distance to a point B 
located an optical distance of unity from point A will be smaller 
than the distance D-E between two points with an average den- 
sity between them, because the path A-C has a higher density. 
This short path will affect all photons coming from outside the 
dense region, giving a linear correction in pqkqcL. Because /(x^, /i) 
are sampled from a smaller volume, the first moment of / which 
corresponds to the flux, will be smaller. 

shows us that a simple formula (given by eq. [ll| ) should de- 
scribe the average density at the sonic point. Moreover, in 
both limits, the unknown constants depend only on geomet- 
rical factors such as the size of the perturbations in units 
of the pressure scale height (/3) or the ratio H between the 
interface surface between high and low density regions and 
their volume in units of d, as well as on the particular shape 
the inhomogeneities adopt. 

In both limits, the unknown constants A~^ and B were 
shown to be proportional to Ar/fcopod — I3~^ . Since the 
perturbations are expected to be of the order of the scale 
height, this factor is of order unity as well. If the pertur- 
bations are smaller, the term will be larger and with it the 
density. 

We define now a wind function VV(r) such that: 

Peonic = ^W(r)punit(r - 1) (28) 

(the factor ^/u is introduced so that the ensuing wind mass- 
loss relation will not include it). From our analysis thus far, 
we found that if the geometrical properties of the developed 
inhomogeneities are not a function of T, then the prefactor 
to (F— 1) should be constant (but might be different for the 
two regimes: F ^ 1 and F — 1 ^ 1). 

Since we have no information yet as to if and how the 
geometrical properties change with F, we assume in the re- 
mainder of the paper that they are constant and that the 
wind function is a constant as well. Thus, 



W(F) = -^{A~^ or B) ^ const. 



(29) 



We should emphasize that a simplifying assumption 
that we implicitly made is that the same proportionality 
relation exists for the large and the small optical depths. 
This is found to be the case in the Markovian mixture of 
low and high density regions. However, this need not be the 
general case. 

We should also bare in mind that the transition be- 
tween the two normalizations (i.e., a different W for the 
optically thick and thin cases) should take place at F — 1 ~ 
1 — Fcff,oo Si 1 — Fcrit ~ 0.15 — 0.5 (depending on the type 
of instability that governs). This will actually be below the 
typical F's we will obtain in the analysis. That is, we are 
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mainly concerned with the optically thick limit. We should 
also be cautious because the function ?/)(Ar) could in prin- 
ciple also be a function of F (and therefore the constants A 
and B as well) , especially if the geometrical behavior of the 
inhomogeneities changes with F. For example, if the typical 
optical depth of the 'blobs' is F dependent. Only a detailed 
numerical analysis could shed more light on this problem. 

The next step in the analysis presented here is to fol- 
low the nonlinear characteristics of the inhomogeneities. 
This is an intrinsically difficult problem since it requires the 
radiation-hydrodynamical understanding of the saturation 
properties of the instabilities in the medium. This analysis 
has commenced in the form of a detailed hydrodynamic sim- 
ulation. A different approach which will soon be completed, 
is to find empirically W and any possible F dependence. 
This can be performed by solving for the super-E ddington 



counterpart to the core-mass luminosity relation (Paczyn- 
ski, which should describe the steady state of novae 

in their post-maximum decline. Then, by comparing to no- 
vae, both Fcff.oo and W can be extracted for novae with 
different luminosities. 



2.4 The mass loss rate 

Once we have estimated the density at the sonic point, we 
can proceed to calculate the mass loss rate. It is given by: 



'sonic ^sonic, where the relevant speed of sound 
is «sonic = \/vVa- This Speed is generally not the adiabatic 
speed of sound since the flow isn't adiabatic. 
Using eq. (]2^), we find: 

. r,2 47rGAf,W(r)(r- 1) 

M — 47r_Ko/5sonicWsonic — 



KqVs 



>V(F)(L-/:Edd) 



(30) 



Note that £Edd becomes the modified Eddington lu- 
minosity if the underlying opacity is not that of Thomson 
scattering. Eq. ( |3c| ) is the basic result of the present theory. 
From the analysis of the location of the critical point in 
we have seen that W should be of order unity or some- 



2.3 



what larger if the size of the perturbations d is similar to the 
pressure scale height of the atmosphere (since W scales with 
Ip/d). Moreover, since W depends only on the geometrical 
properties of the perturbations, it could be close to a con- 
stant if these properties are not a function of F. Hence, we 
assume for simplicity that VV(F) is a constant of order unity. 
We iterate once more that it is only with a more elaborate 
simulation or with more accurate and careful observations, 
that a more accurate functional form can be deduced. For 
the meantime, we will have to settle with this simplifying 
yet reasonable assumption. 

We will show in the rest of the paper that eq. ( po| ) pro- 
vides an explanation to the mass loss from bright classical 
novae as well as from r) Carinae and allows us to connect be- 
tween the observed luminosity and mass loss rate, a relation 
which hitherto did not exist for super-Eddington systems. 



2.5 High Load winds 

Unlike normal stellar winds, the thick winds formed in the 
super-Eddington flows not only have a high mass momen- 



tum relative to the total radiative momentum, which is ex- 
pected in any thick wind, but more importantly, the kinetic 
energy flux can be comparable to the radiative luminosity. 
It also inadvertently implies that the rate of gravitational 
energy 'pumped' into the outflowing matter is also compa- 
rable. Consequently, the luminosity L in eq. (^o|) is not the 
observed luminosity at infinity. Instead, we have to substi- 
tute Lo for L where 



it, 



+ 



(31) 



(The approximation neglects the kinetic energy at the sonic 
points, i.e., assuming that <^ wfac which gives Ltot ~ ^o)- 
The effect of a lowered observed l umino sity was coined 'pho- 
ton tiring' by Owocki & Gayley (1997), who solved for the 
behaviour and evolution of the wind. Their solution related 
the variables Loo, Lq, M and v'^^^. The basic equations 
are the equations of continuity, momentum conservation and 
energy conservations: 



M 
2llr 



const = Anr vp 
GA/*(1-F(r)) 



L{r) = Lo 



M 



(32) 
(33) 

(34) 



where F(r) is the ratio between the luminosity L{r) and 
the local modified Eddington luminosity jCEdd which we as- 
sumed to be constant, since the opacity is assumed to be 
given by the Thomson scattering. The solution to this set of 
equations, af ter neglecting the speed of sound at the base of 
the wind, is (Owocki & Gayley, 1997): 



m 



exp(— ma;)) , 



(35) 



where x is defined as 1 — Ro/r. Ro is the radius of the sonic 
point of the wind and can be described as the 'hydrostatic' 
radius of the star. Below it, the envelope is expanding with 
a subsonic speed, m is defined as MGMt/£,EddRo- It is a 
slightly different definition for the 'photon tiring n umber ' 
than the original definition of m by Owocki & Gayley ( 1997 ) . 
The two are related through rh = Fom. 
At very large radii (r ^ oo), we have 



m 



- exp(— m)) . 



(36) 



Using the wind model developed here (eq. we can 
close the relation between M (or rh) and the luminosity 
(Fo/^Edd) of the star: 



2 Ct, 



-(Fo-l) = (Fo-l)W^, 



(37) 



where W is the 'scaled' wind parameter. This is a dimen- 
sionless version of eq. ( po| ) - the basic mass-loss luminosity 
relation of SEWs. 

For given Fq, Z^Edd and Ro (or v^^^), we can now cal- 
culate rh (or M), and from it calculate Woo (or «^) and 
Foo. To get the latter, we substitute the energy conservation 
equation into the equation for w to get: 



Fo exp(— m) 



(38) 
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In order to compare with observations, we need to translate 
Loo, M and Vod, into an observed temperature. This is done 
using a steady state optically thick wind model. When in- 
tegrated upwards, a photosphere should be obtained where 
the optical depth is of the order of 2/3[|. The general case is 
far from trivial because the opacity is ~ ^KscKab where Kab 
is the absorptive opacity which is generally much smaller 
than the scattering opacity. The simplest approximation is 
first to assume that the opacity is that of Thomson scatter- 
ing. This was done when the original steady s tate optically 



Sha- 



thick wind explanation to novae was proposed ( Bath 
viv. It yields an effective temperature that satisfies 

the following equation: 



log 



M 



gr s-i 



= 22.0 -2.0 log 



+ log 



108 



\W °k) 



0.5 log 



(39) 



IQiS g]-g g-1 



In reality though, the absorption opacity is lower and the 
photosphere is deeper in the wind, where the temperature is 
higher. Thus, a temperature estimate based on eq. (^) will 
be adequate only when we use an effective temperature that 
is defined through aT^g = L/A-kt^y,, where rph is the radius 
where the optical depth for the continuum becomes 2/3. 
This temperat ure is used, for example, in the analysis by 
Schwarz et al. (1998). If a real colour temperature is needed 
(such as by comparison to a Planckian spectrum), then a 



more extended analysis, as was carried out by Bath (1978) 
is needed. Since we do not require a very accurate mass 
loss temperature relation (as the measurement error in the 
observations is much larger anyway), a simple fit to Bath's 
results is sufficient. One obtains that for Tcoiour ^ 6000°K, 
the relation is: 



log 



M 



= 22.3-1. 



, 1 / ^colour \ 

'^"HToJ^J 



(40) 



+ log( 



108 



Tj+0-651og( 



It provides an estimate to logM that is better than 0.1, 
which is more than sufficient for our purpose here. 

Both eqs. (|9|) & ( po| ) can be written more generally as 



M ^ {voo/voo){L^/LY^ 



(41) 



where Ti is either Tofi or Tcoiour and M,Voa,L,T,pL and 
Pt are constants read off eqs. (js^) & (0)- This expression 
allows us to treat both cases simultaneously. 

The last relation that we need in order to close the 
set of equations is the exact value of Vs at the base of the 
wind. To this goal we need the temperature at the base 
of the wind. For simplicity, we assume that the velocity is 
constant and equal to the terminal value; in other words, 
most of the acceleration takes place deep in the optically 
thick part. This assumption is valid as long as rph 3> -Ro- 
under these conditions, the temperature at the bottom can 
be estimated to be: 



t The photosphere in a spherical geometry does not sit at exactly 
T = 2/3, though we assume so for simplicity. 



To 



3/4 



1/4 



(42) 



When we work with the colour temperature, rph is the depth 
of the last mean emission (i.e., t — 2/3 for an opacity of 

\/^^sc^^ab) . 

On one hand, the mass loss depends on the parame- 
ters of the sonic point (e.g., Ro or Vcsc, and To). On the 
other hand however, these parameters depend on the pho- 
tospheric conditions (e.g., through eq. tt2|). Therefore, the 
parameters of the wind, the conditions at the sonic point 
and the conditions at the photosphere should all be solved 
self consistently. It is carried out as follows: 

Assume that the mass of an object M*, its Edding- 
ton parameter at the base To, and a scaled wind parame- 
ter W are given. To obtain the complete wind, sonic point 
and photospheric parameters we begin by calculating m us- 
ing eq. (|37[). Next, Too (and Loo) are found using eq. (Q. 
Woo is then obtained using eq. (bq). To obtain the remain- 
ing parameters, we have to combine together the rest of 
the equations. These are eqs. ( ^) , (ji^), an d (^l|) . The lat- 
ter corresponds to either eq. ^39| ) or eq. (Uoh^depending 
on whether we work with the effective temperature Teff or 
the colour temperature Tcoiour. To close the relations we 
need three more equations, which are Vgsc = 2GM*/i?o, 



ph 



Vs = y^^kTo/^iTTiH = Vs^To/T and crT* = Loo/47rr 
When combined, we find the following expressions for T, 
Woo and the mass loss M: 



Ti 



A-jv^a^G^ L^PL rh^M- 



Edd^ 



i oo i-Edd ^" 



3 + 5p^ 



3 + 5p2^ 



M 



'2C^+^l^^mMwV^£^ 



STrsG'^'-C^dd"^ Af^T''WV2 



c8«fw8rorj 



3 + 5p2^ 



(43) 



The expressions for the additional parameters can be found 
in a straightforward manner. 

The contour plots of roo(ro,W^), i'oo(ro,M^) and 
Tcoiour (ro, W) are depicted in Fig. ^ for a IMq system. 



3 APPLICATION OF THE SEW THEORY 

We now proceed to apply the SEW theory to systems that 
clearly exhibited super-Eddington outflows over long dura- 
tions relative to the dynamical time scales. 



3.1 Two super-Eddington nova 

Although many novae display super-Eddington luminosities 
for at least short periods, it is often difficult to unequivocally 
show that a particular nova was indeed super-Eddington for 
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Figure 4. The state solution in the - W plane. Here, Fq is 
the Eddington parameter at the base of the wind, while W is the 
scaled wind parameter describing the relative 'load' of the wind 
W = Wv^g^/2vsC. Larger W's imply a smaller base radius. The 
thin continuous set of lines is lines of constant Too > 1- The sepa- 
ration between two consecutive lines is 0.2 dex. The set of dotted 
lines is lines of constant Too < 1 separated by 0.2. The Too = 1 
line separating between the two sets of lines, is marked heavier. 
The dashed set of lines is lines of constant Voo separated by 100 
km s"'^. The asymptotic line of ti^o = 0, beyond which there is 
no steady state wind solution, is marked as a thick line. The dot- 
dash set of lines describe constant colour temperatures separated 
by 2500K. Note that the lower lines are the hotter ones. The evo- 
lution of FH Ser is marked by stars connected with a continuous 
line. The evolution of FH Ser from top left to bottom right was 
calculated by combining the following equations as explained in 
the text: (a) The wind model to relate M to Fq (given in this 
work), (b) The temperature mass loss relation (given by Bath 
1978), and (c) Relations between the observed parameters and 
the parameters at the base of the wind (given by Owocki & Gay- 
ley 1997). The stars are the values obtained at different days from 
the data of Friedjung (1987). The fact that the light curve con- 
spicuously follows an iso-tioo contour arises because the observed 
Voo for FH Serpentis happens to vary only little, from 670 to 770 
km s^-"^. 



a long duration, even if it wa^. The reason is that during the 
'bolometric plateau' that a nova exhibits, the luminosity can 
be close to the Eddington limit (from either above or below 
it) such that even small uncertainties in distance and lumi- 
nosity can hide the true nature of the flow. Moreover, since 
the wavelength of maximum emission shifts to the UV dur- 
ing the 'bolometric plateau', as the temperature increases, it 
is difficult to acquire accurate bolometric luminosities from 
Earth based observations. Our insistence of using good data 
results with having only two clear cases in which long du- 
ration super-Eddington flows were observed and have good 
enough data. The two cases are Nova LMC 1988 #1 and 
Nova FH Serpentis. Other cases such as V1500 Cygni are not 
as clear cut, even though one would suspect that long du- 
ration SEWs could have been present. In V1500 Cygni, the 

t On average, the absolute magnitude Mb of a classical novae 
of which th e WD mass 0.5Mq, peaks at a super-Eddington 
luminosity ( |Livio, 1992 ). 



nova was observed to shine at significantly super-Eddington 
luminosities for a few days after the eruption, and when ob- 
served again on day 100, it was 2/3 of Eddington. Thus, 
there is no clear indication for the duration of the super- 
Eddington episode, only a lower limit. Another recent ex- 
ample is Nova LMC 1991, which was observed to be super- 
Eddington for two weeks (Schwarz et al., 2000), reaching a 



truly impressive luminosity as high as 2.6 ± 0.3 x 10 erg 
s"M As we shaU soon see, both V1500 Cygni and LMC 1991 
were classified as very fast novae and hence are marginally 
useful candidates for a steady state analysis. Another nova, 
V1974 Cygni 1992, had very beautiful bolometric observa- 
tions (Shore et al., 1994) which could have been potentially 
very useful for the analysis done here . Unfortunately, th e 
distance uncertainty of 1.8 to 3.2 kpc (Paresce et al., 1991:) 
corresponds to more than a factor of 3 uncertainly in the 
bolometric luminosity, which is too large. 



3.1.1 Nova LMC 1988 #1 

Nova LMC 1988 #1 occurred, as the name suggests, in 
the LMC. Thus, its distance is known more accurately 
than the distance of galactic novae. This distance, coupled 
to an ext ensive multi-wavelength campaign by Schwarz et 
al. (1995), resulted with their rather accurate finding that 
Nova LMC 1988 #1 had an average bolometric luminosity 
of (3.0 ± 0.3) X 10^* erg s"^ during the first 45 days after 
visual maximum. By considering that the maximum mass a 
WD can have is 1.4 Mq, Schwarz et al. (1998) concluded 



that the nova had to be super-Eddington for the long du- 
ration. We use their data and results for the analysis. This 
includes the evolution of the bolometric luminosity and the 
effective temperature (as opposed to a colour temperature). 



3.1.2 Nova FH Serpentis 

Nova FH Serpentis is less clear than Nova LMC 1988 #1. 
Specifically, an accurate distance determination was ob- 
tained onl y after using the HST measurement of the expand- 
ing ejecta (Gill & O'Brien, 2000). The reddening which is in- 
accurately k nown (and is probably in the range E{B — V) = 
0.64 ± 0.16, iella Valle et al., 1997) then remains as the 



main source of error in the bolometric luminosity determi- 
nation. If we take the determination of th e poss ible range 
for the bolo metric luminosity by Duerbeck ( 1992 ) and della 
Valle et al. (1997) and correct for th e som ewhat larger dis- 
tance determined by Gill & O'Brien (2000), we find that the 



bolometric luminosity on day 6.4, when the first UV mea- 
surement was taken, is Mboi = — 7™5±0™2, namely, L must 
have been more than 2.5 x lO"'* erg s~^. 

Furthermore, since the total mass ejected has a lower 
bound of 5 X 10~^ -^0i ^-^d is probably more of the order 
of 2 X lO""* Mq, we should take the kinetic energy of the 
outfiow into account. Since the ej ecta is moving at about 
500 km s" ' ' 



(GiU & O'Brien, 



200C) 



the ejecta has a mini- 
mum kinetic energy of 2.5 x lO** erg. A lower estimate to 
the contribution of the kinetic energy to the total fiux at 
the base of the wind would be to divide the kinetic energy 
uniformly over the duration of the bolometric plateau that 
lasted about 45 days. This gives a minimum contribution of 
ikin,min = 6.5 X 10'^^ erg s"'^. Siucc the total mass loss is 
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probably higher than 5 x 10~^ Mq and since the mass loss 
rate was likely to be higher than average when the luminos- 
ity was higher than the average luminosity, it is likely that 
Lkin > Lkin.min On day 6.4. When we add the kinetic energy 
to the observed bolometric luminosity at day 6.4, we find 
that Ltot at the base of the wind must have been at least 
3.1 X 10^* erg s'^ Clearly, even if the mass of the WD is 
large (which is inconsistent with the eruption not being a 
very fast nova), the luminosities are super- Eddington. The 
fact that the kinetic energy necessarily implies that the nova 
was super-Eddington (at t he ba se of the wind) was already 
pointed out by Friedjung ( 1987 ) who did not have observa- 
tions on the mass ejected but instead used photospheric con- 
straints on M. Here, we have reached the same conclusion 
without resorting to the photospheric analysis. This could 
be done due to the better distance measurements (which 
are within the error but on the high side of previous esti- 
mat es), a nd better reddening measurements of della Valle et 
al. ( [1997| ). 



We use the data of Friedjung (1987, also in Friedjung ^ 
198s|)^r the behaviour of the temperature, luminosity and 
velocity at infinity. These data are b ased on the UV measure- 
ments by Gallagher & Code (1974) and in addition include 
the integrated flux that falls outside the UV range observed 
by Gallagher & Code assuming a black body distribution|^ 
We correct the luminosities to in clude the better distance 
determination by Gill & O'brian (|200(]|). 



3.2 Is a steady state model appropriate for the 
'bolometric plateau' of novae? 

A steady state model is probably appropriate for two rea- 
sons: 



(i) Circumstantial evidence: Kato & Hachisu ( 199^) com- 
pared their steady state model f or win ds to the results of a 
dynamical evolution of Prialnik (1986) and found relatively 
good agreement, implying that a steady state solution is 
valid for most of the nova evolution. 

(ii) Physical reasons: A wind with a terminal velocity of 
Voo originates from a typical radius ro that has an escape ve- 
locity of order Voo ■ Two conditions should therefore be satis- 
fied for the wind to be in a steady state. First, the sub-sonic 
region beneath the base of the wind should be acoustically 
connected on time scales shorter than the typical evolution 
time of the system. Namely: 



^sub 



TO 



< to 



(44) 



The second criterion is that the time it takes the acceler- 
ated material to reach the photosphere is shorter than the 
time at question: 



^sup 



cr — sonic ■ 



^Gvolutio 



(45) 



where tevoiution is the typical time scale for changes 
in the system. Both Nova FH Serpentis and Nova 



§ Note that when taking this extra flux outside the observed UV 
range, the slow increase followed by a decreas e in b olometric lu- 
minosity, as described by Gaflagher & Code (1974), turn into a 
slow decrease of the bolometric luminosity. 



LMC 1988 #1 satisfy both requirements from days 4 and 
2 onward respectively. Thus, a steady state solution for the 
objects ought to be found for tcvoiution ~ month. As previ- 
ously mentioned, Novae V1500 Cygni and LMC 1991 do not 
satisfy the required conditions - they are very fast novae 
that evolve dramatically on a time scale of days, so a steady 
state wind cannot be assumed in their analysis. 



3.3 Application to Nova FH Serpentis 

Using the described metho d of analysis, we proceed to fit 
the data of Friedjung (1987) which includes Tcoiour, Lao and 



Voo- Given a wind parameter W, a white dwarf mass Mwd 
and a hydrogen mass fraction X, we solve for, at each ob- 
served point, the values of Fq and W which result with the 
observed Voo and Too- This is done iteratively. That is, Fq 
and W are assumed to be known such that Too and Uoo can 
be calculated according to the procedure described at the 
end of §2.1:. Then, Fq and W are varied until the calculated 
Foo and Voo agree with observations. (A function of Fq and 
W is defined such that its value is the error between the ob- 
servations and the obtained values of Foo and Voo - The func- 
tion is then minimized using a maximum 'downhill' gradient 
method). 

The resulting set of Fq and W for each observation point 
is then used to predict a colour temperature. The predicted 
temperatures are then compared with the observed temper- 
atures (and their error) to give a number. 

The results for X = 0.3, 0.7 are depicted in Fig. |. We 
find that the possible 1 — cr range for W is about 21.5 to 
34. If we allow larger statistical variation of 3 — cr then the 
value of W can range from 20 to about 39.5 (if we restrict 
ourselves to a reasonable WD mass range between 0.6 and 
1.0 solar masses, as the decay rate would suggest it is). 

An additional systematic error arises from the inaccu- 
racy of the derived reddening of Nova FH Serpentis. This 
uncertainty is portrayed by the error which could increase 
or decrease the derived W by roughly 4-1.8 or -7.2 respec- 
tively. 

Fig. P depicts the observed colour temperature evolu- 
tion and the predicted colour temperature evolution using 
the best fit values of W, Foo(t) and W{t), and the nominal 
values of X = 0.7 and Mwd = 0.8Mq. Good agreement be- 
tween the predicted temperature using the SEW theory and 
the observed colour temperature can be obtained, clearly 
demonstrating that the theory is consistent. 

The temporal evolution of the observed luminosity, as 
given in the same figure, can be integrated to obtain a total 
mass loss of 8 x lO""" Mq. The main source of uncertainty 
in this figure is again the reddening which could increase 
or decrease the mass loss by about -1-10% and -30% respec- 
tively. This result is completely consi stent with the obs e rva- 
tional range of ( 2 - 20) x 10"^Mq ( [Hack et ah, 1993 
h O'Brien, 20001). 



GiU 



An important point which should be considered is that 
the wind is likely to be clumpy, a fact which the analysis thus 
far did not take into account. In the region around the sonic 
point, the optical depth of typical clumps is expected to be 
of order unity. Above this point, typical clumps are there- 
fore expected to be optically thin. Thus, the main effect of 
the dumpiness is to increase the effective absorptive opac- 
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Figure 5. The fit of the predicted temperature (from the lu- 
minosity and wind theory) to the observed temperature variations 
of Nova FH Serpentis, for two values of the Hydrogen mass frac- 
tion X. The 1 and 3 a statistical variation contours are marked 
in the W - MwD plane. Good fits are obtained for reasonable 
parameters. A of 4.2 is obtained for 8 — 2 = 6 degrees of free- 
dom. The systematic errors are large and arise from the inaccurate 
knowledge of Mwd, X and the reddening to Nova FH Serpentis. 
The limit on the systematic variation of the latter is given by the 
arrows. 
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Figure 6. The observed colour temperature behaviour Tobs of 
Nova FH Serpentis (filled triangles, taken from Friedjung 1987) as 
compared with the predicted colour temperature behaviour Tpr 
from the observed luminosity and wind model (open triangles, 
slightly offsetted to the right), for the best fit model assuming 
MwD = O.8M0 and X = 0.70. The additional plots are of the 
mass loss (the stars and dashed line). The mass loss integrates to 
8 X 10~^A/q. Finally we plot as filled squares the ratio of Too to 
Fq . Apparently, about 40% of the original radiative luminosity at 
the base of the wind is used-up under these parameters to acceler- 
ate the wind and compensate for its potential energy, dumpiness 
will reduce this fraction. 



ity (since one will then obtain that Kab oc \/{f^ > \/ (p)^)- 
The effective opacity which is proportional to the geometric 
mean of the scattering and absorptive opacities will increase 
as -\/C where C is a dumpiness parameter: 



C ; 



(46) 



Moreover, from Bath & Shaviv ( 1976| ) we have that for a 
fixed temperature, outflow velocity and luminosity, the in- 
ferred mass loss rate is proportional to k~^^^. Thus, a first 
guess would be that by assuming homogeneity we over esti- 
mate W by roug hlyafactorofC^/*. A more detailed analysis 
which repeats the fitting actually shows that the power law 
relation is over estimated by about 10-15%. That is to say, W 
is over estimated by about C° '^'' for Mwd ~ (0.8 - 1.0)Mq. 

We do not know what the value of C should be. A rough 
guesstimate would be to take the value found in another type 
of system in which strong optically thick winds are observed, 
namely, WR stars. In binary systems in which one of the 
stars is a WR, several measurements of the mass loss in the 
winds can be performed. Some (which are usually the more 
difficult ones) measure the actual mass loss M (for example, 
using scattering induced polarizations or measurement of 
the period slowdown of the binary, P) while others measure 
MC (for example, using radio measurements of the free-free 
emission). In the case of V444 Cygni, a value of C « 3 is 
obtained ( St. -Louis et al., 1993| ). If we adopt this value as a 
typical one, then the inferred range oiC'^^^W, which is 12.5 
- 41.5, corresponds to W = 5.5 — 18. 



3.4 Application to Nova LMC 1988 #1 

Nova LMC 1988 #1 has better data than Nova FH Ser- 
pentis in several respects. First, the absolute luminosity is 
known with higher accur acy. S econd, the temperature ob- 
tained by Schwarz et al. (1998) is an effective temperature 
and not the colour temperature. Consequently, the analysis 
does not depend on the absorptive opacity and therefore the 
dumpiness of the wind. The main draw back however, is that 
Nova LMC 1988 #1 does not have adequate measurements 
of the evolution of the velocity of the outfiow. The veloc ity 
of 1800 ± 200 km s" 



adopted by Schwarz et al. (1998) is 
based on the measurements of optical emission lines that are 
also consistent with the subsequent observations of 'Orion' 
emission lines when the nova was optically thin. There are 
no available records to velocities derived from the principal 
(absorptive) spectrum. 

An important question should be raised. Is the veloc- 
ity of the the photospheric material given by the velocity of 
the diffuse enhanced spectrum or by the principal spectrum? 
The answer is probably the latter. Irrespective of the theo- 
retical argumentation for which spectrum is formed closer to 
the photosphere, the circumstantial evidence points to the 
principal spectrum. This is because the velocities obtained 
from the principal spectrum are similar to those obtained 
by measurements of the nebular expansion years after the 
eruption. Namely, the principal spectrum provides the ve- 
locity of the bulk of the material while the diffuse enhanced 
spectrum probably gives the velocity of the small amount 
of material that is initially injected at high speeds. In Nova 
FH Serpentis for example, the principal spectrum gives ve- 
locities in the range 670 to 770 km s~^. The diffuse enhanced 
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Figure 7. The visual and bolometric luminosity of Nova 
LMC 1988 #1 taken from Schwarz et al. (1998). Open circles 
are the visual luminosities. The measured bolometric luminosi- 
ties are marked by x's and an error bar. Because of the moder- 
ate decay rate of the nova, its mass should be of order or less 
than about 1.2Mq. X = 0.3 is the typical hydrogen mass frac- 
tion in nova ejecta. The dashed line is a least squares fit to an 
exponential decay for the bolometric luminosity, giving Lbol = 
9.3 X lO"*!/© exp(-V90.5 days). The solid line is the best fit for 
a luminosity behaviour for an exponential decay plus a constant, 
giving: Lbol = 4.13 X IO^Lq exp(-t/20.0 days) -1-5.64 X 10*Lq. 
Evidently, the nova was super-Eddington for a long duration. Un- 
fortunately, adequate measurements of the evolution of the veloc- 
ity as a function of time do not exist. 



spectrum is in the range 1300 to 1900 km s-\ Clearly, the 
principal spectrum is closer to the observ ed expansion ve- 
locity of the nebula at 490 ± 20 km s"^ ( |GiU fc O'Brien 
200(|~] 



In order to obtain an estimate for what the principal ve- 
locity is, we tak e the two relat ions that relate tz to ^principal 
and ^diffuse (e.g., Warner, 198!-) and relate the two velocities. 
The 1800 ±200 km s"^ observed for the emission lines, then 
corresponds to a principal spectrum of 900 ± 150 km s~^. 

Taking the above considerations into account, the anal- 
ysis proceeds in a manner similar to that of Nova FH Ser- 
pentis with a few notable exceptions. 

In addition to W, Mwd and X, we leave the velocity 
Hoo as a free parameter and let it vary from 750 km s~^ to 
2000 km s~^, corresponding to the complete range of possi- 
ble velocities, without actually any prejudice for or against 
which spectral velocities should be taken. We assume the 
velocity to be constant with time, which towards the end of 
the eruption, when small mass outflows are expected (and 
therefore a smaller observational footprint), can be a bad 
approximation. We therefore fit the data only to the first 
month of observations. 

Since the dimensionality of the free parameter space is 
too large, we do not plot the contours of but instead find 
and plot the value of W which minimizes it given Mwd , X 
and Woo- Moreover, since the observed luminosity exhibits 
rather large variations, as is apparent in Fig. ^, a smoothed 
functional behaviour is adopted. A fit of the form Lboi = 
aexp (— 6t) + c yields the best fit. 



Figure 8. The best fit values of W for Nova LMC 1988 #1 found 
for different values of the H mass fraction X and different assumed 
constant velocity for the outflow, as a function of the assumed WD 
mass. Values ranging between 5 and 15 are possible. 



The result of the aforementioned procedure are depicted 
in Fig. H Without any prejudice for Mwd, X or Voo (except 
for limiting them to a reasonable range of 0.6Mq < AfwD < 
1.4M0 , 750 km s"^ < Woo < 2000 km s'^ and 0.3 < X < 
0.7), one finds that the wind model can adequately explain 
the results with a wind parameter in the range: 



Wlmc = 10.0 ± 5. 



(47) 



This result is consistent with both the one obtained from 
Nova FH Serpentis and the one expected from the SEW 
theory for (3 somewhat smaller than unity and/or ct either 
closer to 1 or much closer to than the nominal value of 
a = 1/2. 



3.5 The great eruption of rj Carinae 

The case of 77 Carinae is quite different from the novae ana- 
lyzed. First, the systems are entirely different. Instead of a 
solar mass type WD, i] Car is a blue super-giant with a mass 
of order of 60 to IOOMq . While the super-Eddington phase 
of novae lasts of order a few months during which there is 
mass loss with a typical rate of 10"'' Mq yr~^, the giant 
eruption of rj Car which started around 1840, lasted for 20 
years and exhibited a mass loss rate of order 10~^ Mq yr~^. 

A second difference appears in the way the analysis pro- 
ceeds. While the novae analyzed had a detailed evolution of 
the temperature and luminosity, all that we have for r) Car is 
an estimated average luminosity during the 20 year eruption 
and an estimate for the integrated mass loss. Yet we apply 
the theory and stress that the very same theory is applicable 
to a wide range of systems. 



Following Davidson (1999), we adopt the following pa- 



rameters for the eruption of r] Car: Duration of 20 years. 
Integrated radiated flux of lo'*^-^*" '' erg. Ejected mass of 
1 to 3AfQ, and a mass velocity at infinity of 650 km . 
The mass adopted is 80 ± 20Mq which can correspond to 
the estimate of ~ 60Mq if it is a part of a double star and 
~ IOOMq if it is single. 
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To obtain W from the observables, we first assume 
that the mass loss rate was constant and equal to M — 
(1 ± 0.5)Mq yr~^. Given the aforementioned nominal val- 
ues for the luminosity, velocity at infinity, stellar mass and 
mass loss rate, we can iterate for the values of To and W 
according to the procedure outlined at the end of §2.5 and 



the fitting do ne t o each point of Nova FH Serpentis in the 
beginning of §3.3. The difference between the cases of rj Car 



and Nova FH Serpentis is that here we only have one data 
point (corresponding to the 'averaged' eruption) and that it 
does not include a temperature measurement. Therefore, we 
do not have enough observables to obtain a St- Namely, 
we have the correct number free of parameters to fit the 
data exactly. The set of nominal values for the observables 
therefore results with a particular value for W. To obtain the 
error in the estimate, we take each observable and change it 
to the extremes allowed by its error bar. This gives the error 
in W induced by that particular observable. The errors from 
the different observables are then added in quadrature. 
The wind parameter that we find is 



W = 16.2 ± 12. 



(48) 



The greatest contribution to the error arises from the in- 
accuracy of the average luminosity while the second largest 
contribution, which is half the size of the first, comes from 
the inaccuracy of the mass ejected. 

Although the errors in W are larger than those obtained 
for Nova LMC 1988 #1, the range obtained from rj Car in- 
cludes entirely the range obtained from Nova LMC 1988 #1. 
Namely, the results are consistent and the value obtained 
from Nova LMC 1988 #1 should be taken as the best esti- 
mate for W. 



3.6 Possibility of Super-Eddington fluxes 

Although we cannot say at the moment anything quantita- 
tive about the steady state luminosity that will be attained 
without proper knowledge of the properties of the porous at- 
mosphere, we can understand why super-Eddington fiuxes 
are a natural result. 

The steady state luminosity at the post-maximum of 
novae is generally supposed to be given or at least approx- 
imated by the core-mass luminosity relation (CMLR, e.g., 
Tuchman & Truran, 1998) which describes systems in which 
a burning shell is situated on an inert (and hence fixed) core 
(Paczynski, 1970). The CMLR however, saturates at the Ed- 
dington luminosity. So, how can it describe cases which were 
clearly super-Eddington? 

A clear physical understanding of the CMLR can be 
found in Tuchman et al. ( 1983 ) who studied a system com- 



posed of an inert core on top of which there is a burning shell, 
a radiative layer with sharp gradients and a convective layer 
on top of that. They showed that the conditions above the 
radiative layer are unimportant for the determination of the 
burning luminosity (thus, the mass of the envelope is unim- 
portant). Only the burning layer and the radiative layer on 
top of it affect the luminosity. How would a porous atmo- 
sphere change that? Since any inhomogeneities are expected 
to form in the top part of the atmosphere (above the con- 
vective lawr) , it does not change the analysis of Tuchman 
et al. (1983). That is the case only as long as a consistent 
solution for the top part of the atmosphere is obtainable. 



However, as the mass of the WD increases and with it 
the Eddington parameter, the top part of the atmosphere 
becomes unstable against formation of inhomogeneities and 
once these are formed, allows a larger radiative flux for 
the same average temperature and density gradients. If this 
larger flux is close enough (or larger than) the Eddington 
luminosity, one o f the main assumptions in the analysis of 
Tuchman et al. (1983) breaks down because the radiative 



layer ab ove the burning shell necessarily has to become con- 
vective ( Joss et al., 1973 ). Thus, an additional branch for the 
mass-core luminosity relation becomes possible in which the 
top part of the atmosphere is inhomogeneous and below it, 
all the way down to the burning shell, the envelope is convec- 
tive. Since the luminosity now depends on the photospheric 
conditions (the convective layer adjusts itself to 'relate' the 
conditions below and on top of it), the core-mass luminosity 
relation can lose its insensitivity to the parameters of the 
envelope. 

3.7 The 'Transition Phase' 

One of the seemingly odd behaviour displayed by a large ma- 
jority of the classical nova eruptions is a transition phase. 
If it appears, it starts once the visual magnitude has de- 
cayed by 3 to 4 magnitudes. During the transition phase, 
the light curve can display strong deepening, quasi-periodic 
oscillations, erratic changes or other complicated behaviour. 

The origin of the transition phase is not clear and more 
than one explanation had been suggested. For example, the 
transition phase roughly corresponds to the stage when the 
photosphere has shrunk to the size of the binary separation 
such that the companion star can stir up the envelope to 
produce non trivial behaviour. 

The SEW theory naturally introduces another explana- 
tion for the transition phase. If we look at the Fq - W tra- 
jectory of Nova FH Serpentis in Fig. ^, one cannot avoid the 
extrapolation of the trajectory into the zone of 'no steady 
state configuration'. What is this region? 

The mass loss rate is determined by the luminosity. The 
mechanism that generates this wind does not depend on 
whether the luminosity is sufficient or not to carry the ma- 
terial from the sonic point to infinity. If the radius of the 
sonic point is too small and the escape velocity too high, 
the wind simply stagnates before lea ving t he potential well. 
As predicted by Owocki & Gayley (1997), no steady state 



solution for the wind exists under such conditions. In other 
words, the wind model predicts that no steady state could 
be reached if the radius falls too quickly. This appears to 
be a likely scenario in the case of Nova FH Serpentis if one 
extrapolates its trajectory seen in Fig. ^, and indeed. Nova 
FH Serpentis was observed to have a transition phase in 
which the luminosity faded dramatically. 

Both the fading and the erratic or quasi-periodic be- 
haviour seem like a natural result here. From Fig. ^, it is 
apparent that Too can fade dramatically before entering the 
transition phase. This arises from the fact that close to the 
transition phase almost all of the fiux is used to pay for the 
potential toll. Once in the 'domain of no steady state' the 
non trivial 2D or 3D flows that must result could potentially 
result with non trivial variability. 

The faded and/or variable luminosity phase is expected 
to end when the luminosity at the base of the wind falls 
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below the Eddington limit, shutting off the SEW, at which 
point the 'naked' white dwarf should emerge. 



3.8 General Mass Loss of Novae 

We have previously treated two specific novae for which the 
temporal evolution of the luminosity and temperature is 
known in detail. We now turn to describe the nova pop- 
ulation in general for which only the average behaviour is 
known. To do so, we create a template nova as a function of 
white dwarf mass and explore its properties. Its predicted 
mass loss is then compared with the observed integrated 
mass loss. Clearly, we expect the theoretical prediction to 
provide the guide line to which the average observed be- 
haviour should be compared to. 

Let us first estimate the average mass loss in novae. 
To this goal, we use the average trends in the functional 
behaviour of the luminosity, decay time and velocity as a 
function of white dwarf mass to predict the integrated mass 
loss during the super- Eddington episode. Specifically, we use 
the following: 

(i) We use the average relation between t^, the time it 
takes for the visual luminosity to decay by 3 ma gnitudes, to 
the mass of the WD, as is given by Livio (1992). 



(ii) We assume that the photospheric velocity is given by 
the velocity of the principal spectrum. This velocity has an 
average rela tion t o the decay time ta, as given for example 
by Warner ( |l989[ ). 

(iii) We take the peak bolometric luminosity to be the 
peak visual luminosity. This is permissible since novae's 
maxima are generally in the visual. M„ is relat ed to tz, and 
therefore to Mwd , by relations given by Livio ( |l992| ) . Using 
A/wD, Vao and Too, we find the Eddington parameter Fq at 
the base of the wind during the pe ak b rightness. We do so 



using the relations given in section 2.5, 

(iv) We assume that the bolometric luminosity at the base 
decays exponentially (or the magnitude linearly) . That is, it 
has the form: L{t) = Fo-CEdd exp(— f/r). Since we expect the 
transition phase to arise when the luminosity approaches the 
Eddington luminosity, and since the transition phase usually 
sets in after the visual decayed by 3 to 4 magnitudes, or 
about 3.5 on average, we can relate to the exponential 
decay constant of the base luminosity r: 



(3.5/3)t3 
In(Fo) ■ 



(49) 



Few words of explanation are in order. It is conceivable that 
as the nova eruption proceeds, since the amount of fuel is 
fixed and all of it burns simultaneously, we do not get a 
strict steady state cigar-type burning and fixed bolometric 
luminosity. The bolometric luminosity must change grad- 
ually with time. It is customary to assume as a first ap - 
proximation that Lboi = constant ( Bath fc Shaviv, 1976[ ). 
However, a more realistic treatment that more accurately 
describes the observations, is to assume a gradual, though 
slow, decline expressed as an exponential decay with a long 
time scale (see for instance Fig. |7|). For comparison, we also 
repeat the whole calculation using a linear decline of the 
form L{t) = roi2Bdd(l — i/i") and show that the exact form 
of decay is not critical. 

(v) If we integrate the mass loss rate given by the wind 
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Figure 9. Nova ejected shell mass vs. mass of WD. The symbols 
describe the measured mass loss for different observed novae as 
compiled by Hack et al. (1993). Large measurement errors are 
apparent in cases where more than one measurement exists (in 
which case a vertical line connects the points). One should also 
note that some measurements overestimate the mass loss when 
the wind is clumpy, as it is predicted to be. The solid line describes 
the mass loss obtained by the wind model assuming a nominal 
value of W = 10.0 as obtained from Nova LMC 1988 #1, following 
the procedure given in the text. The long dashed line is obtained 
when changing the Hydrogen mass fraction from 0.7 to 0.3. The 
upper (and lower) short dashed lines are obtained when taking a 
value for W which is higher (lower) by 5.0, while the dash-dotted 
line arises when the assumed exponential decay is replaced by 
a linear decay. The dotted lines arise when taking into account 
the natural scatter in the - AfyvD relation obtained by Livio 
(1992). 



theory (eq. |3C|) , we find that the total ejected material during 
the super-Eddington episode is: 



Me- 



Mdt : 



CVs 



)dt 



W 

/lEddT (To - 1 - In(Fo)) , 



(50) 



where Ie '^s the time it takes the bolometric luminosity to 
decline to the Eddington limit. 

The results are plotted in Fig. ^ as a function of the 
mass of the white dwarf, together with the observed de- 
terminations of ejecta masses, as was complied by Hack 
et al. (1993). We first notice that in those cases where multi- 



ple determinations were obtained, there are large differences 
between measurements. These differences should therefore 
be taken as the typical 'error' in cases where only one de- 
termination was performed. The different theoretical pre- 
dictions correspond to changing the wind parameter W to 
within the possible range of 10.0±5, when assuming a linear 
instead of an exponential decay, and when taking into ac- 
count the natural scatter in the - Mwd relation of Livio 



( [l992|) . 



To within the large uncertainties in the ejecta mass de- 
terminations (which have a 1 — cr scatter of 0.54 dex), the 
prediction and observations appear to be in good agreement. 
The super-Eddington episode of novae could account for the 
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bulk, if not all, of the ejected material as a function of Mwd- 
One should nevertheless note that the logarithmic average of 
the observations: ( log ^ q (Mojoct a /Mq yr"^)) = -4.40±0.09 
is on the upper side but still within the allowed region pre- 
dicted by the SEW theory (which has only a small functional 
dependence on A4"wd)- Since some of the ejecta mass deter- 
minations are susceptible to dumpiness, taking the latter 
into account should tend to reduce the logarithmic average. 
If C « 3 is a typical value and if characteristically about a 
half of the measurements actually measure A/ejoctaC, then 
the average of the mass loss would be smaller by about a 
0.25 dex. 

Several additional conclusions can be drawn from Fig. ^ 
The plateau in the predicted mass loss, where the mass loss 
does not depend on the mass of the white dwarf, extends 
from A/wD = O.SAf© to 1.2A^0. This appears to agree with 
observations. Moreover, the extrapolation to WD masses 
beyond 1.2Af0 agrees with the mass loss observed from 
the very fast nova V1500 Cygni, however, this extrapola- 
tion should be done cautiously because the assumption of a 
steady state wind is not strictly valid. 

The plateau predicted by the SEW and seen in the ob- 
servations is counter to theoretical predictions of the TNR 
process in which the general trend should be a smaller mass 
loss with larger WD masses. This theoretical trend arises 
from the fact that more massive WDs are more compact and 
ignite the TNR after significantly less material is accreted. 
For some reason, this trend is not seen. The current TNR 
theory of novae tends to predict values which are about half 
an order to an order of magnitude smalle r than the typi 
cal observations for large WD masses (e.g., St arrfield, 1999 
Prialnik fc Kovetz, 1995| |Starrfield et al., 1998| ). The imple 



mentation of the SEW theory in immorical simulations of 
thermonuclear runaways in now underway with the purpose 
of finding the luminosity at the base of the wind self consis- 
tently as well as to see to what extent the incorporation of 
the SEW changes the predicted amount of ejected material. 

3.9 A 'Constant Bolometric Flux' 

Let us return to Fig. ^ Recall that lines of constant radius 
are close to being vertical. One can see that for moderately 
'loaded' winds and To of a few, an evolution in which the 
temperature increases but the apparent luminosity remains 
constant is possible if the radius does not decrease dramat- 
ically. Under such conditions, the luminosity at the base of 
the wind does fall off to get a higher and higher temper- 
ature. However, the lower mass loss predicted implies that 
less energy is needed to accelerate the material to infinity 
and so a larger fraction of the base luminosity remains after 
the wind has been accelerated. 

This could explain for example the constant bolomet- 



ric lu minosity observed for V1974 Cygni 1992 ( Shore et al, 
1994|) and su pports the early working hypothesis ( [Bath & 



Shaviv, 1974 [Gallagher fc Starrfield, 197^ ) that novae evolve 



with a bolometric luminosity which is almost constant, or at 
least one that does not vary dramatically. 

Interestingly, depending on the load of the wind, the 
constant apparent bolometric fliix can be either super or 
sub-Eddington. That is, even if an apparent sub-Eddington 
luminosity is observed, the system could still have been 
super-Eddington over dynamically long durations. In such 
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Figure 10. Comparison of the wind parameter obtained from 
three independent objects. In the case of FH Ser, the result 
C0.67yy constitutes an upper limit since any dumpiness in the 
wind will mimic a larger inferred value for W if homogeneity 
(i.e., C = 1) is assumed. The value of C = 3 is a guesstimate from 
the dumpiness observed in WR winds. Since the wind is expected 
to be clumpy, the larger value obtained is a good indication that 
the analysis and the model are consistent. The vertical dotted 
lines describe the theoretical prediction and its uncertainties us- 
ing the simple 'blob' model. For nominal values for the unknown 
geometrical parameters (a = 1/2. /3 = 1 and H = 2), the cen- 
tral line is obtained using eqs. @ and (^. The two additional 
dotted lines represent the uncertainties arriving from the geomet- 
rical parameter H describing the surface to volume ratio in units 
of the size of the blobs, which could reasonable vary from H 1 
to S ~ 3. The arrows represent the change in the predicted W 
when changing the volume fraction of the dense regions. When 
large volume asymmetry is present, the value of W increases. It 
is minimized for a = 1/3. /3 can shift W in both directions. The 
dashed vertical line represents the prediction using the Marko- 
vian results (eq. ^). Clearly, the results from the three objects 
are consistent with each other. They are somewhat larger than the 
theoretical prediction using nominal values (/3 = 1, a = 1/2) but 
they do agree with the predictions using geometrical parameters 
well within reasonable bounds (e.g., a ~ 0.8 or larger). 



cases, the kinetic and potential energies of the flow are im- 
portant players in the total energy budget. 



3.10 dumpiness of the wind 

dumpiness was already mentioned on several occasions. 
Since it is important, it deserves a more concentrated dis- 
cussion. 

SEWs are expected to be always clumpy as they are 
generated by an inhomogeneous atmosphere. The results for 
C^'^^yV for Nova FH Serpentis as compared to the results 
for W of Nova LMC 1988 #1 do indicate that the winds are 
clumpy and that the dumpiness factor could be similar to 
that already observed in the optically thick winds of Wolf- 
Rayet stars, namely, C ~ 3. This can be seen in Fig. ^ which 
depicts the W obtained in the three discussed systems. 

dumpiness is also important because it tends to offset 
the estimates for the mass ejected from novae and other ob- 
jects by overestimating it. This arises because many emission 



© 2001 RAS, MNRAS 000, 



Super Eddington Winds 17 



processes are more efficient at higher densities. For example, 
the op acity per unit mass of free-free emission is propor- 



tional tap. Thus, if the material is clumpy, then the higher 
density material is more efficient in producing the observed 
radiation than the case of the material spread evenly over 
the entire volume. 

What is the expected size of the clumps? Since they 
originate from the inhomogeneities at the base of the wind, 
one needs to know the typical size d of the atmospheric 
perturbations at this point. For the instabi lities found t o 
operate in Thomson scattering atmosphere ( Bhaviv, 2001 ), 
the typical size is of the order of the size of the pressure scale 
height Ip (say, pip). As the wind expands from its base, the 
clumps will keep their angular extent relative to the star. 
Thus, the spherical harmonic £ at which the structure should 
peek should be of order: 



Adding t o eg. (po[) the basic resu lts for optically thick winds 
eg. TO (jOwocki fc Gayley, 19971) and eg. Q ( |Bath fc Sha" ' 



2n 



R 

pip 



7r(l - rcff.oo) v^, 



/3~0(1) (51) 



where is the speed of sound at the base of the wind, 
which is higher than that at the photosphere. The typical 
number obtained for novae, is ^ ~ 100. This is a large num- 
ber and it implies that the ejecta has many small clumps 
in it. Since the perturbations are expected to be dynamic 
and change over a sound crossing time, their vertical ex- 
tent should be of order {vao/vs)l3lp, at large distances, after 
the velocity dispersion had time to disperse the blobs ver- 
tically. The horizontal dispersion arising from velocities of 
order of the speed of sound are not important because of the 
horizontal expansion induced by rarefaction in the spherical 
geometry. 



4 DISCUSSION & SUMMARY 

The existence of steady state super-Eddington outflows is 
an observational fact. One therefore needs to explain on one 
hand how atmospheres can sustain a super-Eddington state 
and on the other, one needs to understand the winds that 
they generate. 

We have tried to present the following coherent picture: 
Homogeneous atmospheres becomes inhomogeneous as the 
the radiative flux approaches the Eddington limit. This is 
due to a plethora of instabilities. The particular governing 
instability depends on the details of the atmosphere. As a 
conseguence of the inhomogeneity, the effective opacity is 
reduced as it is easier for the radiation to escape, and con- 
seguently, the effective Eddington limit increases. 

Super-Eddington configurations are now possible be- 
cause the bulk of the atmosphere is effectively sub- 
Eddington so that the average is also sub-Eddington. Very 
deep layers advect the excess total luminosity above Ed- 
dington by convection. Higher in the atmosphere, where 
convection is inefficient, the Eddington limit is effectively 
increased due to the reduced effective opacity. The top part 
of the atmosphere, where perturbations of order of the scale 
height become optically thin, has however to remain super- 
Eddington. Thus, these layers are pushed off by a continuum 
driven wind. 

By identifying the location of the critical point of the 
outflow, one can obtain a mass-loss luminosity relation. The 
relation, given by eg. (BOf) is the main result of the paper. 



viv, 1976D or eg. (|40D ( path, 1978| ), we are left with one free 
universal dimensionless function W. It was shown that W 
depends primarily on the geometrical characteristics of the 
inhomogeneous atmosphere from which the wind emerges. 
Having no detailed information on the geometrical proper- 
ties, we assumed that they are independent of F, so that W 
becomes a constant. 

To check these ideas, we analyze 2 novae as well as the 
massive star 77 Carinae for which sufficiently detailed and 
accurate data is available. Although the two types of sys- 
tems are notably different, as they have masses, luminosities 
and mass loss rates which differ by orders of magnitude, the 
wind mass loss and the wind parameter are found to be in 
agreement with the theoretical expectation: 



M 



W^^^(F - 1), with W = 10.0 ± 5. 



(52) 



The evolution of the temperature predicted from the lumi- 
nosity agrees well with the temperature measured directly 
in the two novae. 

Another interesting agreement is the consistency with 
clumping. Clumping in SEWs is a natural prediction since 
the atmospheric layers beneath the sonic point are predicted 
to be inhomogeneous. Moreover, dumpiness is a necessary 
ingredient in the present theory that allows super-Eddington 
luminosities. The present theory predicts therefore, that 
SEWs are clumpy. In the analysis of Nova FH Serpentis, 
the wind parameter obtained is coupled to the dumpiness 
of the wind since a clumpy wind with the same observed 
colour temperature will have a lower inferred mass loss. The 
lower M reduces the measured wind parameter towards that 
obtained for Nova LMC 1988 #1 when the typical dumpi- 
ness factors seen in WR winds are taken into account. 

The consistency of all the results is best demonstrated 
in Fig. |lo| which shows that W obtained from the three dif- 
ferent objects is consistent with one another. Although W 
obtained is somewhat larger than the theoretical prediction 
using nominal values for the geometrical characteristics, the 
uncertainties in the latter are large enough to comfortably 
accommodate the range of 'measured' W's. 

An additional way of viewing the results is as follows. 
The SEW theory predicts a range of W which due to un- 
certainties can span at least an order of magnitude. The 
analysis of Nova LMC 1988 #1 yields a W which agrees 
with this range. The analysis of Nova FH Serpentis yields a 
W which too agrees with this range if dumpiness, which is 
predicted, is included. Next, instead of finding the values of 
W which best fit 77 Car, we could have asked the opposite 
guestion. Using the SEW theory which was 'normalized' us- 
ing Nova LMC 1988 #1, we could have calculated what the 
observed mass loss should have been according to the SEW 
theory, and find out that its prediction is completely consis- 
tent with the observed mass loss, even though the object's 
characteristics are different from novae by several of orders 
of magnitude. 

We also identify the occurrence of the 'transition phase' 
observed in more than two thirds of the novae with the ad- 
vance of the atmosphere into the 'no steady state region'. 
As the nova explosion progresses, its luminosity and radius 
decline. However, if the radius decreases too guickly, at some 
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point the predicted SEW will be too heavy for the luminos- 
ity at the base to push mass to infinity. No steady state will 
then exist. The inconsistency might explain the strange be- 
haviours observed in the transition phase of different novae. 

The wind model presented is by no means a complete 
theory for novae since it cannot predict ab initio the lumi- 
nosity at the base of the wind. To obtain the latter, one 
needs to solve for the complete evolution of a nova taking 
into account the fact that the nova's atmosphere becomes 
porous and feels a reduced effective opacity. One expects 
that the lowered opacity increases the luminosity obtained 
in the core-mass luminosity relation, and super-Eddington 
values therefore arise naturally. Preliminary numerical solu- 
tion of the stellar structure equations appear to confirm this 
prediction. 

In order to fully understand the behaviour of the atmo- 
sphere that drives these winds, a 2D or even 3D numerical 
simulations on scales of the order of few optical depths are 
unavoidable because of the intrinsic nonlinear properties of 
the problem. These simulations are underway and will al- 
low one to study in detail the physics controlling the mass 
outflow. 
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APPENDIX A: THE EFFECTIVE OPACITY 
WHEN NOT FULLY IN THE OPTICALLY 
THICK OR THIN LIMITS 

We proceed here to calculate the corrections to the effective 
opacity in systems which are neither fully in the optically 
thick limit nor fully in the optically thin limit. This analy- 
sis is presented to strengthen the heuristic arguments given 
in the section calculating the location of the critical point. 
It also depicts the problems arising when facing the task of 
calculating the radiation transfer in a nonlinear inhomoge- 
neous medium. We begin with an analysis of a system that 
is not fully in the optically thick limit and continue with an 
estimating the effective opacity when not fully in the thin 
limit. 

Al Regions close to the optically thick limit 

Although it is apparently impossible to solve in a close ana- 
lytical form the opacity reduction in an arbitrary nonlinear 
conflguration, we attempt to solve it in at least one realistic 
configuration in order to estimate more accurately B and 
Ph- We shall now assume for simplicity that the perturba- 
tions are vertically elongated. This is actually a reasonable 
approximation because the unstable modes found in Shaviv 
typically have ~ Ip^ but a somewhat smaller kz. 
This implies that most of the interface between high density 
and low density regions is vertical and that the isothermal 
contours are mostly horizontal. 

Since the medium is close to the thick limit, we solve 
for the radiation in the vicinity of the interface, assuming 



(pool) 
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that it is flat and the two density regions extend infinitely 
to either side. 

The framed region shown in figure |^ corresponds to a 
typical area which we wish to solve. We assume that both 
sides of the region have B = Bq + bz. The dense region on 
the right has ph = {1 + S)p and the underdense on the left 
has Pi = (1 — 5)p. 

Using the formal solution given by eq. ( p^ , we can write 
the radiation intensities. For points located in the high den- 
sity region, at a horizontal coordinate xo, the intensity for 
trajectories coming from the positive direction is given by: 



r/ - \ r, , _tan/icos(/p 
I(xo,n-) = Bo + OKop- 



(l + <5)(x- 



cos p 
xo) exp 



(Al) 



Kop{x - Xo)(l + S) 



COS p 



dx, 



where p is the direction of the ray relative to the x-axis 
(perpendicular to the interface) such that s = s/cos/i and 
z = a:: tan ^ cos ^. For trajectories coming from the negative 
x-direction, we have: 



/(x, n+) = Bo + hKop 



tan p cos 



cos p 



(A2) 



(1 + 5){x — Xo) exp 



i^op{x — xo){l + 5) 



+ 



(1 — 5){x — Xo) exp 



cos p 

Kop[2;o(l + 5)+a;(l-(5)] 



cos p 



dx 



dx 



This time cos p is negative. If we now integrate to get the 
flux, we find 



F{xo > 0) = 



47rfe 



+ 



(1 



(A3) 



3cKop(l + 5) ' Scpno 
X [xo5{l + 5) (-6 + io(l + <5)^)r(0,io(l + 5))) 
+ (5(-4 + £o(l + S){-1 + xo + xoS)) exp(-io(l + S))] 

with xo = xoKop. r{n,x) is the incomplete F function (and 
should not be confused with the Eddington parameter F). 
The first term is the uncorrected flux Fh in the high density 
region while the rest of the expression is the correction near 
the interface with the low density region. The net correction 
to the average flux of the system from only corrections in 
the positive x-axis is: 

Kotio 2cKop{l~S){l + 5y^ 



(A4) 



The flux in the negative half, where the density is p(l — 5) 
is obtained by the transformation 5 — + ~5. The sum of the 
net correction to the flux on both sides of a:: = is the total 
correction to the average flux in the system: 

7v(S/Vtot)b5^ 



AF = AF+ + AF- = - 



CKop{l — S^)^ 



(A5) 



The average flux of the system, had it been homogeneous, 
is -Fb = 47r6/3cKopo. If we then assume F S> 1 or 5 — > 1 and 
use eq. (p^), we find: 



AF 



35 



32q(1 



^KOPoVtc 



(A6) 



This result should be compared with the one obtained 
in the more heuristic treatment (eq. The correction to 
the flux obtained here has an additional factor of 3/ (8a) 



but is otherwise the same. That is, if we compare KcH as 
given by eq. (^) with Tp — BAt~'''\ we obtain ph = 1 and 
S = (3H/32Q(l-a)2)//3. 



A2 Regions close to the optically thin limit 

We begin with systems approaching the thin limit. To calcu- 
late the radiation transfer we should perform a linear anal- 
ysis in which the zeroth order terms are those of the thin 
limit. That is to say, the zeroth order terms are those arising 
when the optical depth of the perturbations can be neglected 
(even if 5p/p~ 0(1)). 

The formal solution to the radiation transfer is given by 
eq. dig]). It assumes that the boundaries are many optical 
depths away. This is valid also in the optically thin case 
because even if a scale height is optically thin, the generated 
wind is optically thick unless F — 1 is extremely smaller than 
unity. We also implicitly assume that the photon crossing 
time is significantly shorter than any dynamical time scale 
in the system on which the perturbations can change their 
structure. 

We write the radiation source function of the unper- 
turbed region as -80(1") and the density as po(r), while the 
corresponding perturbations are -Bi(r) and pi(r). 

If we perturb the formal solution of the transfer equa- 
tion (eq. |1E|) we find: 



7o(x,n) 
/i(x,n) 



XoBoe-'^^'ds 

I poo 

XoBie'''°'ds+ / xiBoe'^'^'ds 
Jo 



XoBoe-"""' xids ] ds. 



(A7) 



The small parameter used for the expansion is Ei = 
Jg Xii-^')^^' , namely, the perturbation to the extinction E 
which appears in the exponent. If Ei is much smaller than 
unity, than exp(— i?o — Ei) « exp(— i5o)(l — Ei). Note that 
-El ^ 1 for s's of order the m.f.p of a photon or smaller. 

Next, we assume for simplicity that the unperturbed ra- 
diation source function can be written as: -Bo(x) — Bo-l-x-b. 
After the linearization, it is meaningful to look for harmonic 
solutions to the equations of the form Xi — Xie^^'^ where 
Xi is any perturbed quantity. 

The expression for the first order perturbation to the 
specific intensity can then be obtained rather straightfor- 
wardly if we change the integration order for the double in- 
tegral and choose an integration coordinate system in which 
b aligns with the 2-axis (as opposed to aligning k with z). 
One then finds: 



/i(x,n) 



Xo-Bi ^ 1 (n-b)xi 
ik ■ fi -I- Xo Xo ik • n 4- xo 



(A8) 



The first term describes changes in the source function due 
to local changes in the temperature. The second term de- 
scribes changes in the optical depth due to temperature and 
density perturbations along the line of sight. 

Once the expression for the perturbed radiative flux is 
known, the moments of / can be evaluated directly from the 
perturbed radiation field. Starting with the lowest moment. 
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J is: 



Ji(x) 



1 

_1_ 

Ai 



/i(x, n)di7 
ik ■ n + xo 



(A9) 



+ 



Xo 



(n • b)xidn 
zk • fi + Xo 



— Si - J— ^Aa 
Xo / 3xo 



Xi_ 

Xo 



To obtain the last equality we choose the integration coor- 
dinate system so that k is parallel to its z-axis. In this par- 
ticular choice the denominators of the integrals obtain their 
simplest form. We also defined the following three functions 
for convenience: 

Ai(C) = |tan-^(e) 



A2(C) 
A3(C) 



3 



■tan-^(0 

3[(e2+l)tan-l{«)-5] ■ 

5? 



Ai(C»l) 

A2(C> 1) 
A3(C » 1) 



TT 1 

2 e 

3 

37r 1 

4 



The third function will be used shortly. The unperturbed 
value of J is Jo = Bq. 

The first moment of / is related to the flux (F = 
(47r/c)H) and given by: 



This fact provides the justification for neglecting the tem- 
perature perturbations in the heuristic picture given in § ^ 
Moreover, careful analysis of the definition of the Ai's will 
reveal that the first two terms actually cancel each other. 

The effective o pacity of the medium can now be calcu- 
lated using Shaviv (1998): 



Xeff = 



(Hx) _ Hoxo + (Hixi) + (Hi) xo + Ho (xi) 



Hoxo + (^iXi) 
Ho 



Ho + (Hi) 



o 2 

o Xo 



(z-k)k)f 



(A13) 



where we have used the ^ 3> 1 expansion of A3 and assumed 
without loss of generality that b is in the z direction. We 
point out that had we worked with the opacity per unit mass 
K instead of the extinction x, which is the opacity per unit 
volume, the averaging () would have been a mass weighted 
average and not a volume average. 

It would be interesting to solve the problem exactly in 
the same case which is solved for in the thick limit, namely, 
that of vertical slabs of width d and a density of po(l ± S). 
To do so, we use the Fourier series expansion for the density: 
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47V 
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/i(x, n)ndf2 

XoABidCl 
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where k is a unit vector in the direction of k. The first 
term arises from perturbations to the temperature while the 
second and third from perturbations to the opacity. 

When close to the Eddington limit, the region where 
the optical depth of a scale height is of o r der unity should 
be in t he N AR limit (e.g., [Shaviv, 200l| , [Clatzel fc Kiri- 
akidis, 1993) unless the absorption opacity is extremely 
small (smaller than about (prad/Pgas5)(ws/c) times the ratio 
between the absorption opacity and the scattering opacity). 
Consequently, one can assume that heating or cooling of the 
gas is instantaneous and that B = J. If we impose Ji — Bi 
(sin ce the zeroth order terms satisfy it anyway), we find from 
eq. dAlCj) that: 
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If we write ^ = (fc/xo) then the correction to the flux is: 
k Alio k-b 
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(A12) 



First, we note that the first two terms are proportional to 
^^'^ for large ^s or short wavelengths while the third term 
falls only as such that it would dominate for low ^s. 



p{x) = po 1 + 



45 sin [(2n + l)7vx/d] 
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For Thomson atmospheres we find that xi ~ 1^0 pi- If the 
identity J2n=o'^/^'^'^ + ^ 7C(3)/8 is used (with 
being the (^-function), we obtain: 
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(A15) 



We can also write 5^ — {k,o — i^cS,cx,)/k,o in the case of vertical 
slabs. If we now compare this result to eq. with ?/>(Ar ^ 
1) Ri 1 — Ar^', we obtain: pi = 1 and A ~ 1.28dfi:opo/Ar — 
1.28/3. 

In the more general case of a 7^ 1/2, an analytic solution 
does exist but it involves poly-logarithm functions of a which 
would appear in A, changing it by 0(1). 
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